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1.  INTRODUCTION 

The  purpose  of  this  work  is  to  describe  a straightforward  geometrical  construction 
that  approximates  the  mapping  of  solar  magnetic  fields  through  the  corona  to  inter- 
planetary space.  The  construction  is  based  conceptually  on  the  source-surface  model 
introduced  by  Schatten  et  al.  (1969)  and  by  Altschuler  and  Newkirk  (1969).  The  source 
surface  in  their  models  is  a sphere  of  radius  r^  = 1.6-2. 5 r^  (solar  radii)  on  which  the 
scalar  potential  V (from  which  B = is  derived  for  r < r^)  is  made  constant.  This 
constraint  has  the  effect  of  making  ^radial  at  r = r^.  The  field  outside  the  source  sur- 
face in  their  models  is  generated  from  the  spiral-field  construction  of  Parker  (1958)  by 
imposing  the  continuity  of  B at  r = r . 

The  traditional  source-surface  models  (Schatten  et  al.,  1969;  Altschuler  and 
Newkirk,  1969)  are  computationally  very  simple  and  yield  a qualitatively  pleasing 
f ' reproduction  of  major  heliomagnetic  features  such  as  helmet  streamers,  filaments, 

r arcades  of  closed  field  lines,  and  coronal  holes.  However,  these  models  are  only 

I approximations  to  a rigorous  magnetohydrodynamic  (MHD)  solution  of  the  interaction 

I between  the  solar  wind  and  the  coronal  magnetic  field.  The  inaccuracy  inherent  in  these 

t models  prevents  a detailed  comparison  of  magnetic  structures  in  the  outer  corona  emd 

I 

solar  wind  with  observed  magnetic  field  patterns  in  the  photosphere.  The  major 
' weakness  in  these  models  is  the  incompatibility  of  the  postulated  radial  magnetic  field 

[ with  MHD  at  r < 10  r^.  This  weakness  is  well  illustrated  by  comparison  with  the  only 

I existing  MHD  solution  (Pneuman  and  Kopp,  1971a,b),  which  was  obtained  for  the  special 

i case  in  which  the  normal  component  of  the  photospheric  field  is  that  of  a dipole  field 
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with  a polar  field  strength  of  1 G.  In  '^’igure  1 the  usual  source-surface  solution  for  this 
case  is  compared  with  the  MHD  solu  ion.  The  differences  are  quite  obvious.  The  MHD 
field  lines  are  not  even  approximately  radial  for  r > r^  = 2.5  r^,  and  the  choice  of  a 
smaller  source  surface  would  make  the  disagreement  even  worse. 

There  have  been  subsequent  attempts  to  replace  the  source-surface  models  by 
models  which  represent  the  MHD  effects  more  accurately.  Schatten  (1971)  developed  an 
analytical  technique  that  was  based  on  the  idea  of  placing  sheet  currents  between 
oppositely  directed  open  field  lines.  The  resulting  construction  produced  a qualitative 
improvement,  in  that  the  open  field  lines  were  indeed  deflected  equatorward  instead  of 
being  radial;  however,  field  lines  in  the  region  near  the  cusps  of  helmet  streamers  were 
poorly  modelled  thereby.  (In  order  to  model  well  the  direction  of  the  open  field  lines,  the 
cusps  had  to  be  placed  unrealistically  low,  i.e.,  at  1.6  r^  in  the  dipole  case.)  Yeh  and 
Pneuman  (1977)  developed  a much  more  rigorous  sheet-current  model  in  which  the 
location  and  strengths  of  current  sheets  between  oppositely  directed  open  field  lines  (and 
between  closed  and  open  structures)  are  actually  calculated  iteratively  by  determining 
the  force  balance  between  the  magnetic  field  and  the  gas  pressure  in  the  expanding 
corona.  Their  result  was  shown  to  be  very  accurate  when  compared  with  the  MHD  dipole 
solution,  but  the  calculation  was  so  lengthy  that  it  appears  impractical  to  generalize 
their  scheme  to  the  point  of  applying  it  to  a set  of  observed  photospheric  magnetic 
fields. 

In  the  present  work  we  develop  a construction  that  retains  most  of  the  compu- 
tational simplicity  of  the  earlier  source-surface  models  but  which  significantly  improves 
the  agreement  with  the  MHD  solution  in  the  dipole  case.  We  do  this  by  adopting  an 
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Fig.  1.  Configuration  of  Magnetic  Field  Lines  Emanating  from 

an  Internal  Solar  Dipole  in  the  MHD  Solution  of  Pneuman 
and  Kopp  (1971a,  b;  dashed  curves)  and  in  the  Source- 
Surface  Model  of  Altschuler  and  Newkirk  (1969;  solid 
curves).  The  spherical  source  surface  (dotted  curve) 
was  taken  to  have  a radius  equal  to  2.5  rg 
(Altschuler  and  Newkirk,  1969). 


appropriate  non-spherical  source  surface  on  which  the  helio magnetic  scalar  potential  V is 
made  constant.  The  exterior  solution  is  constructed  in  this  case  by  mapping  the  outward 
normal  to  our  source  surface  in  accordance  with  a generalized  Parker  spiral.  Our  choice 
of  a non-spherical  source  surface  is  well  motivated  by  MHD  considerations,  and  the  field 
that  we  generate  agrees  well  in  magnitude  and  direction  with  the  one  calculated  by 
Pneuman  and  Kopp  (1971a, b). 
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2.  NONSPHERICAL  SOURCE  SURFACE 


For  our  non-spherical  source  surface,  we  select  an  appropriate  isogauss  of  the 
spherical-harmonic  extrapolation  of  the  magnetic  field  derived  from  internal  sources, 
i.e.,  without  considering  the  effects  of  currents  external  to  the  sun.  The  identification 
of  the  source  surface  with  an  isogauss  rather  than  a sphere  represents  a better  simulation 
of  magnetohydrohynamic  (MHD)  behavior.  One  should  expect  the  magnetic  field  to  be 
distorted  very  little  by  plasma  currents  where  B is  large,  but  very  much  where  B is  small. 
Thus,  the  influence  of  the  solar  dipole  should  be  relatively  less  perturbed  by  plasma 
currents  at  polar  latitudes  than  at  equatorial,  for  observers  stationed  at  the  same 
heliocentric  distance. 

We  do  not  attempt  to  specify  the  source  surface  self-consistently,  e.g.,  ty 
including  the  effects  of  external  currents  on  JB^,  and  therefore  on  the  shape  of  the 
critical  isogauss.  We  are  not  attempting  here  to  solve  a problem  in  MHD  rigorously,  but 
are  instead  trying  to  formulate  a straightforward  geometrical  construction.  The 
accuracy  of  our  construction  is  easily  tested  by  comparing  our  results  with  an  MHD 
solution,  and  our  model  is  amply  vindicated  by  such  a comparison. 

After  the  non-spherical  source  surface  rji  6 , ^)  is  identified  from  the  spherical- 
harmonic  expansion  deduced  from  the  observed  photospheric  magnetic  field,  the  coronal 
field  interior  to  the  source  surface  is  derived  under  the  assumption  that  this  volume  is 
current-free.  The  coronal  B field  there  is  calculated  from  a scalar  potential  of  the  form 
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N n 

^ ‘‘0  2 (gj]™  cos  mv>  + hJJ' sin  m<p  ) PJJ’C  0 ) 

n=l  m=0 


+ n + hJI'  sin  m^  ) Pjj'  (6  )]  . (1) 


where  r^  is  a (constant)  radius  representative  of  the  non-spherical  source  surface.  In  the 
dipole  solution  that  follows,  the  parameter  r^  is  the  equatorial  radius  of  the  source 
surface.  In  a general  solution,  r^  is  merely  a parameter  that  conveniently  represents  the 
particular  isogauss  level  that  is  chosen.  For  example,  one  might  identify  the  critical 
isogauss  either  by  specifying  Bq  (the  value  of  l^lon  the  surface)  or  by  specifying 


fj/r,  ■ ♦ (g\f  * (hJ)2]l/S  (2) 

as  a general  definition  clearly  consistent  with  (and,  indeed,  motivated  by)  the  dipole 
case.  The  coefficients  gj|*  and  hJJ'  are  (of  course)  already  determined  from  the  spherical- 
harmonic  expansion  deduced  from  observations  of  photospheric  field  (e.g.,  Altschuler  et 
al.,  1977). 

The  coefficients  g ^ and  h||'  are  chosen  so  as  to  minimize  the  mean-square 
tangential  component  of^over  the  source  surface.  This  is  an  operational  way  of  saying 
that  B should  be  (as  nearly  as  possible)  normal  to  the  source  surface.  This  statement  is  a 
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straightforward  generalization  from  the  model  characterized  by  a spherical  source 

surface.  There  are  two  main  complications  in  this  generalization:  There  is  no  longer  a 

simple  relation  between  the  g (h^)  coefficients  and  the  g|^  (hJJ^)  coefficients,  and  one 

cannot  be  assured  that  any  spherical-harmonic  potential  exists  that  makes  B normal  to 

the  source  surface  over  the  entire  area.  We  therefore  ask  only  that  the  s ^ and  h'’^ 

^ n n 

coefficients  be  such  as  to  minimize  the  integral 

0-  = y dA  (3) 

over  the  source  surface,  where  B^  is  the  tangential  component  of  B. 

To  obtain  ^ outside  the  non-spherical  source  surface,  we  construct  the  outward 
normal  to  the  source  surface  and  make  the  normal  component  of  B continuous  there. 
This  is  an  obvious  generalization  from  the  case  of  a spherical  source  surface.  It  should 
be  quite  effective  when  the  source  surface  has  a relatively  simple  shape,  e.g.,  in  the 
dipole  case  discussed  in  the  following  section.  However,  one  may  question  this  procedure 
for  more  complicated  source  surfaces,  and  one  must  recognize  that  the  procedure 
becomes  meaningless  if  the  source  surface  is  not  everywhere  convex.  However,  even  in 
such  extreme  cases,  the  solution  internal  to  the  source  surface  might  be  utilized  to 
prescribe  a magnetic  field  geometry  so  that  the  hydrodynamic  problem  of  coronal 
expansion  within  a specified  field  geometry  can  be  solved  (Durney  and  Pneuman,  1975). 
In  such  a case  the  external  flow  field  might  be  obtained  from  hydrodynamics  (rather  than 
MHD)  by  imposing  the  Durney-Pneuman  solution  as  a boundary  condition  at  the  non- 
spherical  source  surface.  In  essence,  then,  we  specify  a very  simple  extrapolation  of  the 
field  normal  to  the  source  surface  but  recognize  that  a more  realistic  (albeit  more 
complicated)  solution  might  be  necessary  in  practice. 
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3.  THE  DIPOLE  CASE 

The  application  of  our  technique  to  the  case  of  an  internal  solar  dipole  is  very 
instructive  for  several  reasons.  First,  the  source  surface  can  be  expressed  analytically, 
so  as  to  make  the  least-squares  minimization  of  (3)  particularly  simple.  This  allows  us  to 
investigate  the  stability  of  the  solution  and  the  convergenge  of  the  calculated  g®  as  N, 
the  maximum  allowed  order,  is  increased.  Second,  the  solution  is  directly  comparable 
with  the  MHD  dipole  solution  of  Pneuman  and  Kopp  (1971a,b).  This  comparison  will  be 
used  as  the  primary  test  of  the  accuracy  of  our  technique.  Finally,  even  such  a simple 
field  configuration  as  the  solar  dipole  already  yields  very  instructive  information  about 
the  behavior  of  the  interplanetary  field,  as  observed  at  r ~ 1 AU.  Therefore,  we 
present  the  detailed  solution  of  the  solar-dipole  case  to  illustrate  the  advantages  of  a 
non-spherical  source  surface.  The  isogauss  surfaces  of  a dipolar  magnetic  field  have  the 
general  form 

r = (1  + 3 cos^e  rg,  (4) 


where  rg  is  the  equatorial  radius  and  6 is  the  colatitude.  For  a 1-G  polar  field  the 
choice  rg  = 2.5  r^  would  correspond  to  a surface  of  constant  B s jjB^|  = 0.032  G in  the 
absence  of  external  currents.  It  develops  (see  below)  that  an  equatorial  radius  Tg  = 2.3  r^ 
provides  a much  better  overall  agreement  between  our  model  and  the  solution  of 
Pneuman  and  Kopp  (1971a,b).  We  remark  parenthetically  that  the  source  surface 
described  by  (4)  encloses  a volume  equal  to  1.380173  (4it  /3)rg,  and  that  a dipolar  source 
surface  of  equatorial  radius  rg  = 2.245  r^  would  enclose  the  same  volume  as  the  spherical 
source  surface  (of  radius  2.5  r^)  used  by  Altschuler  and  Newkirk  (1969). 
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3.  1 Internal  Solution 

In  the  dipole  case,  which  is  azimuthally  symmetric,  the  scalar  potential  inside 
the  source  surface  simplifies  from  (1)  to  the  form 


N 

^ = '‘e  ^1  ® ^ ■'O  8 1 Pn(cos  0) , 

n=l 


(5) 


= 0 


where  the  g ^ are  expansion  coefficients  to  be  determined  below  and  the  Pj,(cos  6)  are 
Legendre  polynomials.  The  term  involving  g®  generates  the  field  of  the  underlying 
internal  dipole.  The  magnetic  field  B {=  -^V)  corresponding  to  (5)  is  given  by 


t _0 


L=  + ^g^Cr^/r^sine 


N 


n(r/r/"^gJp^(cos  0) 


n=l 


N 


0^  ^2  ^ Pj^'(cos  0)  sin  0 . 


(6) 


n=l 


- 0 


The  coefficients  g ^ are  to  be  chosen  so  as  to  minimize  or  , which  is  the  integral  of 
s (fi  X B)^  over  the  surface  defined  by  (4),  where  n is  the  unit  vector  in  the  direction 
of  the  outward  normal  to  the  source  surface.  In  terms  of  (4)  this  means  that  n is  the  unit 
vector  in  the  direction  of  VrJr,  0): 
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Thus,  we  seek  to  minimize  the  integral 


s (£  bJ  dA  s (f  (hx  B)^  dA 


(8) 


with  respect  to  the  f where 


dA  = 2Trrsine  [ r^  + (dr/de  de 


= 2it  r^j  (1  + 7 cos^e  + 8 cos^G  (1  + 3 cos^G  sin  G dG  (9) 


is  the  element  of  area  on  the  source  surface.  Since  the  expansion  coefficients  g ^ in  (5) 
must  vanish  for  even  n in  order  to  preserve  the  north-south  symmetry  of  the  B field,  it 
follows  from  (8)  that 


= 4TTr2 y*  (1  + 


,2.-2/3 


3x^)  (1  + 7x^  + 8x’)  B.'-dx, 


M 1/2 


(10) 


where  x = cos  G . Moreover,  it  follows  from  (6)  and  (7)  that 


= (1  + 7x^  + (1  + 3x^r^^^  (1  - x^)^ 


X I {r^/r/  (1  + x2)  gj 


^ (1  + 3x2)("''2)/6  -0  [nxPj^(x)  + (1  + 3x2)P^'(x)]  \ . 


We  seek  to  minimize  o-  with  respect  to  each  of  the  expansion  coefficients  by 
setting  dr/  9g  ® = 0.  There  follows  from  this  condition  the  requirement  that 


J (1  + 7x^  + (1  + 3x^)'^^  (1  - x^)  (1  + 3x2)('"+2)/6 


X [mxP^(x)  + (1  + 3x‘')P^'(x)]|(r^/rg)^  (1  + x^)  gj 


(1  + 3x2)  ^"^2)/6  + (1  + 3 x2)  P^'(x)]  } dx  = 0 


for  every  odd  integer  m < N.  Thus,  one  must  solve  a set  of  simultaneous  linear  equations 
to  determine  the  required  coefficients  g Our  requirement  can  be  expressed  more 
concisely  in  matrix  form.  We  require  that 
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(13a) 


^nm  = -W^m- 
n=l 


where 


A 


nm 


J (1  + 7x^  + (1  + - x^) 


X [(1  + 3x^)  P^'(x)  + nx  Pj^(x)] 

X [(1  + 3x^)  (x)  + mx  Pji^(x)]  dx  (13b) 


and 


K^=  f (1  + 7x^  + 8xV^^^  (1  + 3x^)  (1  _ x'*) 

A 


X [(1  + 3x^)  P^'(x)  + mxPj^(x)]  dx  (13c) 


for  n,m  = 1,3, 5,7,.  . . , N-2,  N.  The  matrix  dearly  symmetric,  i.e., 

Aj^^  = Therefore,  one  may  alternatively  view  g as  a column  vector  and 

express  (13a)  in  the  form 


N 

V A„„(6S/g?)  = -(r>n)^K 


(14) 


where  is  considered  a column  vector, 
m 

The  least-squares  condition  that  determines  the  g ® can  thus  be  solved  by  standard 
matrix-inversion  techniques.  The  array  is  very  well  conditioned,  being  nearly  tri- 

diagonal. Away  from  the  three  main  diagonals  the  array  values  decrease  by 
approximately  two  orders  of  magnitude  per  diagonal.  We  have  solved  (14)  for  the  g ® 
with  N = 5,7,9,11,13,15.  The  results  are  shown  in  Table  1.  It  is  quite  evident  from  this 
tabulation  that  the  results  are  weU  behaved,  and  that  the  expansion  coefficients  g ^ 
converge  rapidly  to  the  values  that  they  would  acquire  in  the  limit  N — ► . 
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^ ; Table  1.  Optimal  Values  of 

Determined  by  Minimizing  the  Value  of  ^ Specified 

by  (10). 


^ ' \ 

n 

N=5 

N=7 

N=9 

N=ll 

N=13 

N=15 

k 

% 

1 

-0.62917 

-0.63124 

-0.63185 

-0.63205 

-0.63211 

-0.63213 

; "i 

3 

+0.11960 

+0.12537 

+0.12716 

+0.12775 

+0.12795 

+0.12802 

i ' 

5 

-0.03038 

-0.04088 

-0.04448 

-0.04574 

-0.04619 

-0.04635 

t 

7 

+0.01086 

+0.01607 

+0.01820 

+0.01903 

+0.01935 

9 

-0.00418 

-0.00674 

-0.00795 

-0.00847 

11 

+0.00168 

+0.00293 

+0.00360 

' : 

13 

-0.00069 

-0.00130 

15 


+0.00029 


3.  2 Exterior  Solution 


1 


The  meridional  component  of^at  r > (1  + 3 cos^  in  the  present  model 

is  defined  by  constructing  the  outward  normal  to  the  source  surface,  i.e.,  by  constructing 
a straight  line  parallel  to  the  unit  vector  n specified  by  (7).  In  view  of  the  azimuthal 
symmetry  of  the  dipole  problem,  it  proves  convenient  to  introduce  the  cylindrical 
coordinates  p = r sin  6 and  z = r cos  0.  Since  r = p sin  6 + z cos  6 and 
e = p cos  9 - z sin  6,  it  follows  from  (7)  that 

•i 

n = (1  + 7 cos^e  + 8 cos^e  [(1  + 4 cos^G  ) p sin  6 + 4 z cos^el  . (15) 

Thus,  the  outward  normal  to  the  source  surface  at  colatitude  6^  is  a straight  line  having 
the  slope 

dz/dp  = 4(1  + 4 cos^  6 )"^  cos^  6^  ctn  0 (16) 

s s s 

and  the  intercepts 


Pg  = Tq  (1  + 3 cos^  ®s 

(17a) 

Zs  = Tq  (1  + 3 cos^  0g)^^®  cos  0g. 

(17b) 
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3,  3 Comparison  with  MHD  Solution 


In  Figure  2,  the  interior  and  exterior  solutions  are  joined  at  r = (1  + 3 cos^ 
and  are  compared  with  the  MHD  solution.  The  footpoints  of  the  field  lines  were  chosen 
to  coincide  with  those  in  the  figure  published  by  Pneuman  and  Kopp  (197 la, b).  A 

reasonable  match  to  the  overall  configuration  of  field  lines  is  achieved  with  r„  = 2 3 r 

0 ‘ ©’ 

The  improvement  of  the  present  model,  with  its  non-spherical  source  surface,  over 
models  with  spherical  source  surfaces  can  be  seen  by  comparing  Figure  2 with  Figure  1. 
The  direction  of  the  field  lines  in  the  exterior  region  is  improved  dramatically.  Even  the 
shape  of  the  field  lines  interior  to  the  source  surface  is  considerably  improved.  On  the 
basis  of  this  test  with  the  dipole  case,  one  would  expect  the  present  technique  not  only  to 
allow  a much  more  accurate  comparison  of  observed  interplanetary  fields  with 
photospheric  fields,  but  also  to  model  coronal  structures  more  faithfully. 

As  has  been  stated  before,  the  parameter  r^  is  adjustable.  It  is  not  surprising  that 
our  good  overall  fit  to  the  MHD  solution  requires  rg  = 2.3  r^,  whereas  the  equatorial  cusp 
of  the  MHD  solution  occurs  at  r = 2.5  r^.  One  could  not  expect  our  source-surface  model 
to  agree  with  the  MHD  solution  in  every  detail.  Thus,  the  construction  of  a realistic 
coronal-field  model  requires  a certain  degree  of  judgment  in  the  selection  of  r^j,  just  as 
the  selection  of  a radius  for  the  spherical  source  surface  has  required  in  the  past.  Of 
course,  the  parameter  r^  can  be  adjusted  to  match  ..pacific  features  of  the  MHD  solution 
very  well.  Figure  3 illustrates  two  such  adjustments  of  r^:  the  first  to  match  the  MHD 
solution  in  the  region  close  to  the  neutral  sheet,  and  the  second  to  provide  an  explicit 
tracing  of  the  cusped  field  line. 
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Configuration  of  Magnetic  Field  Lines  Emanating 
from  an  Internal  Solar  Dipole  in  the  MHD  Solution 
of  Pneuman  and  Kopp  (1971a,  b;  dashed  curves)  and 
in  Our  Source -Surface  Model  (solid  curves).  Our 


nonspherical  source  surface  (dotted  curve)  is  taken 
to  have  an  equatorial  radius  rQ  = 2.3  r , which 
corresponds  to  a polar  radius  ro  ~ 2.  9 r 

(cf.  Fig.  1).  We  have  taken  N = 15  in  (5). 


Fig.  3.  Configuration  of  Magnetic  Field  Lines  Emanating  from 
an  Internal  Solar  Dipole  in  Our  Source -Surface  Model 
for  Selected  Equatorial  Radii  rQ  > 2.  3 rg,:  (a)  to  Match 
the  Lowest- Latitude  Open  Field  Lines  Shown  by  Pneuman 
and  Kopp  ( 197  la,  b)  at  3 rg  £ r £4  r^,  rp  = 2.  3380  r©; 
and  (b)  to  Illustrate  the  Cusped  Field  Line  Explicitly, 
ro  = 2.3391  r©.  As  in  Figs.  1 and  2,  the  illustrated 
field  lines  from  the  source -surface  model  were 
chosen  so  as  to  intersect  the  sun  (r  = r©)  at  the  same 
set  of  colatitudes  0©  as  the  MHD  field  lines  plotted 
by  Pneuman  and  Kopp  (1971a.  b).  We  have  taken 
N = 15  in  (5). 
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I > 3.  4 Interplanetary  Magnetic  Field 

The  present  model  yields  a calculable  interplanetary  magnetic  field  that  can  be 
compared  with  spacecraft  observations.  This  circumstance  provides  a further  test  of  the 
accuracy  of  the  model  and  allows  for  instructive  predictions  about  the  overall  structure 
of  the  interplanetary  ^field.  Detailed  calculations  are  provided  here  for  the  case  of  an 
internal  dipole. 

Figure  4 illustrates  the  calculation  that  is  required.  We  must  calculate,  as  a 
function  of  heliocentric  distance  r,  the  cross-sectional  area  of  the  magnetic  flux  tube 
that  is  bounded  in  latitude  by  a pair  of  infinitesimally  separated  normals  to  the  source 
surface.  The  flux  tube  in  question  is  a figure  of  revolution  about  the  z axis.  By  virtue  of 
its  construction,  the  flux  tube  contains  a fixed  number  of  field  lines.  Thus,  the  product 
of  (a)  the  cross-sectional  area  of  the  flux  tube,  and  (b)  the  meridional  component  of  the 
magnetic  field,  must  remain  independent  of  r throughout  the  flux  tube. 

The  element  of  area  transverse  to  the  flux  tube  is  given  by  dA  = 2Tr  p dSj^  , where 

ds^  is  the  thickness  of  the  flux  tube  at  heliocentric  distance  r = p esc  0 . However,  the 

thickness  dSj^  at  (r,  e,  is  proportional  to  the  distance  of  the  point  (r,  9,  (f)  from  the 

center  of  curvature  (C)  of  the  source  surface  for  the  point  (r  , e , <p  ) along  the  same 

s s s 

normal  (see  Figure  4).  Thus,  if  R is  the  radius  of  curvature  at  (r  , 0 , ^ ),  the  value  of 

s s s 

ds  is  proportional  to  R + 6 , where  6 is  distance  from  (r  , 0 , ) to  (r,  0,  (f). 

* s s s 

Therefore,  the  meridional  component  of  B at  (r,  Q,  v)  is  related  to  the  normal 
component  B„  of  B at  (r„,  0„,  by  the  formula 
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Fig.  4.  Geometry  of  Meridioneil  Cross  Section 
of  Magnetic  Flux  Tube  (bounded  by 
solid  lines  that  intersect  at  C),  as 
Employed  in  the  Calculation  of  B at 
Points  (r,  P)  Outside  Our  Nonspher- 
ical  Source  Surface  (dashed  curve) 
From  the  Normal  Component 
(Bs=  n • B ) of  B at  the  Source  Surface 
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1 


Bm  = (Ps/p)[R/(R+  6)]Bg 

= (Pg/r)  (sin  e^/sin  6)  [l  + (6  /R)]  "^(1  + 3 cos^  (18) 

We  calculate  R as  a function  of  and  obtain  6 and  6 as  functions  of  r and  6 . It 

s s 

would  have  been  more  convenient  to  obtain  6^,  R,  and  6 in  terms  of  r and  6,  but  this 
proved  intractable. 

Using  the  cylindrical  coordinates  p = r sin  e„  and  = r^  cos  0 , we  find  that  the 

3 s s s s s 

source  surface  defined  by  (4)  satisfies  the  equation 

(Ps  + = (P^  + 4z2)r®  . (19) 


We  consider  z^  to  be  a function  of  p^  on  the  source  surface  and  differentiate  (19) 
implicitly  to  obtain 

,1  2 3 

Zg  = - j (1  + 4 COS  0g)  sec  0g  sin  0^  (20a) 


z''=  - (1/8  z ) [ (1  + 3 cos^  0 J sec®  0 sin^  0 
^ S s s s 


The  radius  of  curvature  is  thereupon  given  (e.g.,  Thomas,  1960)  by 


R=  iz:r^[l  + . (21) 

s s 

We  note  in  particular  that  R-*-oo  at  the  equator  (8^=  tt/2)  and  that  R = 0.8  r^  at  the 
poles. 

Since  z^  is  the  slope  of  the  source  surface,  the  equation 

p - Pg  = -{z-Zg)Zg'  (22) 

characterizes  the  normal  to  the  source  surface.  Therefore,  the  normal  distance  6 of  the 
point  (r,  0 , <p)  from  the  source  surface  is  given  by 


6 - [(P  - p/  * fe  - 


[1  + ( z -Zg)  sgn  (cos  6g).  (23) 


The  remaining  task  is  to  calculate  z - z.  in  terms  of  r and  0 . To  do  this  we  write 

s s 

r^  = p^  + z^  = (p  - Pg+  Pg)^  + (z-Zg  + Zg)^ 

= (p  - Pg)^  + (21-Zg)^  + rg  +2pg(p  - Pg)  + 2Zg(z-Zg)  (24) 
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and  eliminate  (p  - p ) with  the  aid  of  (22).  The  result  is  a quadratic  equation  in  (z  - z ) 
that  yields  two  solutions,  one  of  which  can  be  identified  as  extraneous  because  it  fails  to 
reduce  to  z - z^  = (r  - r^)  sgn  (cos  e^)  at  cos  0^  = ± 1 (^  , at  the  poles).  The  accept- 
able solution  reads 


4[l  + (zp^J  (z  - Zg)  = 

{ 16  [l  + (zp^](r^  - r^  + r^  (1  + 3 cos^  0 )^  sec®  0 sgn  (cos  0 ) 
o 9 s s s s 

-r^  (1  + 3 eos^  0 J sec^  0 . (25) 

O d S 

where  r^  = (1  + 3 cos^  ’’o  usual.  The  substitution  of  (25)  in  (23)  yields  6 as  a 

function  of  r and  0^  for  fixed  r^.  Moreover,  it  follows  from  (23)  that 

cos  0 = (r^/r)  (1  + 3 cos^  9 cos  0 

0 s s 

+ [l  + (zp^]“^^^  ( 6/r)  sgn  (cos  0^)  , (26) 

i.e.,  that  0 can  be  written  explicitly  as  a function  of  r and  0 for  fixed  r„. 

— r so 


I 
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We  identify  in  (18)  with  the  normal  component  of  at  the  source 

surface: 

= (1  + 7 cos^  0.  + 8 cos^  [(1  + 3 cos^  9 J B^  + sin  0^  cos  9^  B_  ] . (27) 

s s s srssv 

In  order  to  facilitate  applications  of  the  present  model  by  the  interested  reader,  we 

3 0 

present  (in  Table  2)  listings  of  the  quantity  (r„/r„)  (B/g. ) sec  0^  as  a function  of  cos 

U 0 S 1 s s 

at  the  source  surface  for  selected  values  of  N in  (5).  We  have  selected  this  particular 

combination  of  variables  because  of  the  fact  that  its  limit  remains  finite  (and 

independent  of  r^/r^)  as  one  approaches  the  equator  (0  = Tr/2),  i.e.,  because 

®s  * ®1  ®s  limit.  The  expression  for  B^^  given  by  (18)  becomes 

indeterminate  at  cos  0 = ± 1,  i.e.,  at  the  poles.  Since  we  find  = 0 and  (therefore) 

R - ^ “ 0.8  r at  cos  6 = ±1,  the  correct  limit  for  (18)  is  obtained  by  considering  a 

s s s 

spherical  expansion  from  the  center  of  curvature,  which  is  located  on  the  axis  of 
symmetry  at  z = 0.2  r sgn  (cos  0).  Thus,  we  obtain 


®m  = ^r  = 16[5(r/r  ) - 1 ] B (28) 


at  cos  0 = ±1.  The  source-surface  radius  at  the  poles  is  given  by  r^  = 2^^^^  rjj, 
according  to  (4).  Since  B^  /B^  = - z^'  along  the  entire  field  line,  it  is  easy  enough  to 
reconstruct  B in  cylindrical  ( p , ^ , z)  or  spherical  coordinates  outside  the  source  surface: 
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Table  2.  Values  of  (rp/r^)^  (S.B/gj)  sec  0^  at  the  Nonspherical  Source 


Surface  Given  by  (4),  as  Computed  from  (5)  with  the 

—0  * 

Optimal  Coefficients  g Specified  in  Table  1 


COS  0^ 

N=5 

N=7 

N=9 

N=ll 

CO 

II 

N=15 

0.0 

3.452 

3.745 

3.960 

4.113 

4.222 

4.300 

0.1 

3.396 

3.674 

3.873 

4.009 

4.101 

4.164 

0.2 

3.243 

3.479 

3.629 

3.717 

3.766 

3.792 

0.3 

3.028 

3.193 

3.269 

3.293 

3.291 

3.281 

0.4 

2.773 

2.844 

2.840 

2.809 

2.781 

2.765 

0.5 

2.484 

2.459 

2.401 

2.362 

2.352 

2.358 

0.6 

2.171 

2.082 

2.032 

2.034 

2.053 

2.063 

0.7 

1.853 

1.773 

1.787 

1.819 

1.823 

1.810 

0.8 

1.572 

1.585 

1.634 

1.624 

1.607 

1.611 

0.9 

1.400 

1.489 

1.455 

1.444 

1.463 

1.455 

1.0 

1.442 

1.245 

1.376 

1.289 

1.346 

1.312 

B s n*B  at  the  source  surface, 
s Mm 


r 

r 

I 


1 


Bp  = - [l  + (zp^  Zg  sgn  (cos  Gg)  (29a) 

= [l  + (zp^  ] Bj^  sgn  (cos  8^)  (29b) 

B = [l  + (z'S^  ] (cos  0 - z ' sin  0 ) B sgn  (cos  0 ) (29c) 

Id  d 111  d 

B = - [l  + (z')^]"^'^^  (sin  0 + z ' cos  0)  B sgn  (cos  0 ) (29d) 

V d d m d 


We  assume  that  the  sun  rotates  with  angular  velocity  about  the  dipole  (z)  axis  and 
that  the  solar-wind  velocity  u is  directed  along  the  normal  to  the  source  surface.  Thus, 
we  neglect  the  azimuthal  component  of  u while  requiring  that  u-  /u„  = B.  /B„.  It 

AAV  or  6 r 

follows  from  the  argument  offered  by  Parker  (1958)  that 

B = - (r  sin  6 - r sin  0 ) ( n /u)  B (30) 

Y s s rn 

— fi  —1 

in  our  (more  general)  model.  We  would  take  n = 2.8925  X 10  sec  in  order  to  repro- 
duce the  27-day  periodicity  of  solar  magnetic  features  observed  from  the  orbiting  earth 
(cf.  Dessler,  1967).  A simple  model  for  B^  could  be  constructed  by  adopting  the 
nominal  quiet-time  solar-wind  velocity  (u  = 320  km/sec)  recommended  by  Hundhausen 
(1970),  but  a better  model  v/ould  require  us  instead  to  specify  the  probable  variation  of  u 
with  r and  0^.  Rather  than  burden  the  present  work  with  a model  for  the  solar-wind 
velocity,  however,  we  shall  ignore  interplanetary  comparisons  of  B^  with  that  obtained 

from  our  model  and  concentrate  instead  of  B_. 

m 

i 

1 
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Figure  5 illustrates  the  variation  of  with  r/r^.  In  the  left  panel  one  follows  a 

field  line  by  holding  fixed  the  value  of  0 (colatitude  at  which  the  field  line  intersects 

source  surface)  or,  equivalently,  of  0^  (colatitude  at  which  the  field  line  intersects  the 

surface  of  the  sun);  in  this  case  the  latitude  of  evaluation  decreases  with  increasing  r/r  . 

0 

In  the  right  panel  one  evaluates  at  fixed  latitude  \ = |(Tr/2)-0i.  Thus,  the  left 

panel  illustrates  a conceptual  point,  while  the  right  panel  illustrates  an  observational 

one.  Both  panels  illustrate  a pronounced  departure  of  (r/r^)  B^  from  constancy  at  low 

latitudes,  in  marked  contrast  to  the  case  of  a spherical  source  surface.  With  a spherical 

source  surface  (Schatten  et  al.,  1969;  Altschuler  and  Newkirk,  1969)  one  would  have 

obtained  0 = 0^,  B.  = 0,  and 
s u 

®r  ” ^ ^1  ® 

2 

with  r„  independent  of  0„.  It  would  follow  from  (31)  that  (r/r  ) B , as  well  as  0, 
s s 0 m 

should  remain  constant  along  a field  line.  Moreover,  it  would  seem  from  (31)  that  the 
equatorial  current  sheet  at  0 = •n-/2  is  not  a sharply  defined  feature.  The  azimuthal  ( (p ) 
component  of  the  current  density  would  be  given  by 

J = (c/4Tr ) (p  • V X B = (3c/4it  r ) (r  /r)^  g?  sin  0 (32) 

T AM*  O O 1 

in  the  case  of  a spherical  source  surface.  This  would  represent  a distribution  of  current 
having  a full  width  of  120°  (out  of  a possible  180°)  at  half  maximum. 
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The  results  of  the  present  model  (with  its  non-spherical  source  surface)  stand  in 
sharp  contrast  to  (31)  and  (32).  Figure  6,  for  example,  illustrates  the  variation  of 
with  cos  e at  r = 215  r (ai  AU),  as  obtained  by  us  from  (18).  We  find  B at  rasl  AU  to 
be  nearly  constant  (within  10%  of  the  polar  value)  over  at  least  85%  of  the  unit  sphere. 
Our  meridional  component  of  B,  in  contrast  to  that  described  by  (31),  actually  passes 
through  a broad  off-equatorial  minimum  at  a latitude  ~30°  and  rises  to  more  than  four 
times  the  polar  value  at  a latitude  of  approximately  0.023°;  we  have  magnified  the  scale 
of  the  abscissa  by  a factor  of  100  for  cos  6 < 0.01  in  order  to  show  the  region  of  field 
reversal  clearly  (see  Figure  6).  The  maxima  in  ! B^  ' are  separated  by  ~ 1.2  x 10^  km  at 

O 

a heliocentric  distance  of  1.5  x 10  km  (1  AU).  Such  precise  numerical  results  should 
not  be  taken  literally,  since  source-surface  models  are  intended  only  to  simulate 
dynamical  calculations  and  not  to  supplant  them.  However,  we  view  the  evolution  of 
such  a sharply  defined  interplanetary  current  sheet  as  a major  improvement  over  past 
source-surface  models  that  employed  spherical  source  surfaces. 

A final  comparison  between  our  model  and  that  of  Altschuler  and  Newkirk  (1969) 

is  shown  in  Figure  7.  Interplanetary  field  lines  in  each  case  have  been  traced  from 

colatitude  e at  r = 215  r («!  AU)  to  colatitude  0 at  the  source  surface  (r  = r ) and  to 

0 s s 

colatitude  6^  at  the  surface  of  the  sun  (r  = r ).  We  selected  r = 2.5  r for  the  case  of  a 

0 0 S 0 

2 1/6 

spherical  source  surface  and  r = 2.3  (1  + 3 cos  0 ) r for  the  case  of  a non-spherical 
source  surface,  in  accordance  with  the  examples  analyzed  above.  We  recall  that  0 and 
0g  are  equal  for  the  case  of  a spherical  source  surface,  but  we  note  with  surprise  the 
nearly  linear  relationship  between  0 and  0^  for  the  case  of  a non-spherical  source 
surface.  Moreover,  except  at  cos  0 = 0 and  cos  0 = ±1,  we  find  that  the  field  lines  at 
interplanetary  latitude  \ = |(tt/2)  - 0|  map  to  higher  latitudes  at  r = r^  and  at  r = r^  in 
our  model  than  in  the  model  of  Altschuler  and  Newkirk  (1969). 
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100  COS0  (at  1 AU 


COS  0 I at  1 AU  I 


Variation  of  Meridional  Component  Bm  o; 
Interplanetary  Magnetic  Field  at  Helio- 
centric Distance  r = 215rG®“  1 AU  with 
Colatitude  0 for  Nonspherical  Source 
Surface  with  Equatorial  Radius  rg  = 2.  3 : 
Dashed  curve  and  lower  sccile  apply  for 
cos  0 ^ 10-2;  solid  curve  and  upper  scade 
apply  for  cos  0 s 10-2.  We  have  taken 
N = 15  in  (5).  The  ordinate  measures  Bm 
in  ■Y(ly  = 10 "5  G)  if  one  takes  g9  = 1 G. 


0 AT  1 AU 


Fig.  7.  Results  of  Tracing  Magnetic  Field  Line 
from  Colatitude  6 atr  = 215r0»  1 AU 
in  Interplametary  Space  to  Colatitude  6^ 
at  the  Source  Surface  (r  = rg)  and  Thence 
to  Colatitude  6q  at  the  Surface  of  the  Sun 
(r  = r©).  Dashed  curves  correspond  to 
the  model  of  Altschuler  and  Newkirk 
(1969),  i.  e . , to  a spherical  source  sur- 
face of  Radius  rg  =2.5  r©.  Solid 
curves  correspond  to  the  present  model, 
i.  e . , to  a nonspherical  source  surface 
with  equatorial  radius  rg  = 2.  3 r©.  We 
have  taken  N = 15  in  (5). 


4.  APPLICATION  TO  REALISTIC  B FIELD 


The  obvious  extension  of  our  concept  of  a non-spherical  source  surface  is  to  apply 
it  to  the  spherical-harmonic  expansion  deduced  from  a set  of  observed  photospheric  B- 
field  data.  This  procedure  carries  the  promise  of  allowing  a calculation  of  realistic 
coronal  and  interplanetary  magnetic  fields  that  are  accurate  enough  to  permit  a detailed 
investigation  of  a variety  of  solar-interplanetary  structures,  e.g.,  neutral  sheets,  sectors, 
high-speed  streams,  coronal  holes,  helmet  streamers,  and  so  on.  In  order  to  do  this,  one 
must  develop  a numerical  program  to  minimize  or  as  given  by  (3)  over  an  arbitrary, 
azimuthally  asymmetric  source  surface.  This  is  a rather  large  and  complicated  task,  and 
is  beyond  the  scope  of  the  present  paper.  However,  before  such  a task  is  started,  we  can 
and  should  investigate  what  the  shape  of  a realistic  source  surface  (as  determined  by  the 
procedure  outlined  in  Section  2)  would  be.  Accordingly,  we  have  obtained  from 
R.  H.  Levine  (personal  communication,  1977)  a set  of  spherical-harmonic  expansion 
coefficients  (g^  and  h|^  from  n = 1 to  n = 20,  m = 0 to  m = n)  for  use  in  (1)  for  two  well- 
studied  solar  rotations  (1602  and  1609).  We  have  used  these  coefficients  to  calculate 
several  isogauss  surfaces.  For  purposes  of  illustration  we  have  adopted  the  isogauss 
|b'  = 0.03  G as  the  hypothetical  source  surface.  In  accordance  with  the  present  pre- 

.%w 

scription,  the  coefficients  and  hj|’  in  (1)  have  been  neglected  in  calculating  the 
isogauss. 

The  results  are  shown  in  Figure  8 in  the  form  of  synoptic  contour  plots  of  constant 
r^/r^.  The  actual  shape  of  the  isogauss  surface  for  rotation  1602  can  be  seen  more 
clearly  in  Figure  9,  where  two  meridional  cuts  through  the  surface  are  shown.  The 
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Fig.  8.  Contours  of  Constant  (specified)  t^/tq  on  the 
Isogaussian  "Source"  Surface  |b(  = 0.3  G for 
Rotations  1602  (a)  and  1609  (b).  Based  on  (1) 
with  (gJJ^,  hJJ^)  Supplied  by  R.H.  Levine  (per- 
sonal communication,  1977)  and  (2)  with 
(irT.  Neglected 


particular  longitudes  of  the  cuts  were  chosen  so  as  to  display  the  largest  possible 
variations  of  r^/r^  with  latitude.  In  general,  the  deviation  of  isogauss  surfaces  from 
sphericity  is  very  large,  which  indicates  that  the  application  of  a non-spherical  source 
surface  (rather  than  a spherical  source  surface)  to  coronal  magnetic-field  models  based 
on  real  data  will  produce  significant  changes. 

The  above  solar  rotations  occurred  during  the  Skylab  era  (rotation  1602  in  June 
1973;  rotation  1609  in  December  1973).  As  a result  there  is  a wealth  of  information 
available  for  these  two  rotations,  and  we  can  compare  details  of  the  shapes  of  the  cor- 
responding isogauss  surfaces  with  other  observed  features  ( cf.  various  references  in 
Zirker,  1977).  The  basic  shape  of  the  isogauss  surface  for  rotation  1602  resembles  the 
shape  of  the  isogauss  surface  for  the  dipole  model,  with  obvious  "bulges"  in  the  polar 
regions.  Superimposed  on  this  basic  dipole  shape  are  various  structures  in  the  equatorial 
region  associated  witli  active  regions,  particularly  between  longitudes  of  0°  and  180°. 

This  is  the  same  general  pattern  that  was  inferred  from  observations  of  coronal  holes  and 

interplanetary  streams.  The  isogauss  surface  for  rotation  1609  displays  a considerably 

different  appearance;  its  shape  might  best  be  described  as  that  of  a tilted  dipole.  Again, 

this  basic  shape  was  inferred  from  interplanetary  and  coronal  observations.  (See  i 

Hundhausen,  1977,  for  a more  detailed  discussion  of  the  coronal  and  interplanetary 

phenomenology  during  this  time  period.)  Viewing  each  isogauss  surface  in  more  detail, 

we  note  that  there  is  a one-to-one  correspondence  between  "bulges"  in  the  surface  and  j 

coronal  holes.  For  example,  the  bulge  centered  at  about  10°  longitude  and  20°S  latitude  , 

in  rotation  1602  is  coronal  hole  number  1.  All  of  these  considerations  lead  us  to  conclude 
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that  the  use  of  a non-spherical  source  surface  would  lead  to  a very  significant 
improvement  in  current  coronal  field  models. 


The  failure  of  a source  surface  to  be  everywhere  convex  entails  no  obvious 
impediment  to  the  interior  solution,  i.e.,  to  the  minimization  of  o-  in  (3).  However,  this 
failure  will  definitely  impact  our  construction  of  the  exterior  solution,  since  outward 
normals  from  different  points  on  the  source  surface  will  necessarily  intersect  somewhere 
outside  the  source  surface.  Such  intersections  suggest  the  formation  of  shocks  in  the 
outer  corona  and  perhaps  correspond  to  features  that  are  actually  present  in  space.  If 
confirmed  by  observation  or  by  MHD  theory,  the  presence  of  such  shocks  would  be  of 
interest  in  the  context  of  heating  the  solar  wind  (C.  F.  Kennel,  personal  communication, 
1978).  However,  since  we  do  not  know  the  fundamental  relationship  (if  any)  between 
source-surface  constructions  and  MHD,  it  would  be  premature  to  conclude  that 
indentations  in  our  source  surface  correspond  to  shock-producing  indentations  in  a 
properly  computed  MHD  critical  surface  (C.  F.  Kennel,  personal  communication,  1978). 
The  effect  of  indentations  in  an  isogauss  source  surface  for  the  present  work  is  to 
complicate  the  construction  of  our  exterior  solution,  probably  by  requiring  a hydrodyna- 
mical  treatment  of  solar-wind  expansion  from  the  source  surface.  Such  indentations 
could  perhaps  be  removed  by  smoothing  without  introducing  a major  error  in  the  field 
configuration.  We  do  not  propose  an  arbitrary  smoothing  of  the  source  surface  as  such. 
Rather,  we  contemplate  diminishing  the  spatial  resolution  with  which  ^is  described  at 
r = r^  (e.g.,  by  averaging  the  observed  JB^  over  a larger-than-necessary  area  of  the  solar 
surface)  so  as  to  reduce  the  influence  of  the  higher-order  spherical  harmonics  in  the 
vicinity  of  r = r . Such  a procedure  would  tend  to  preserve  the  gross  shape  of  the  source 

9 

surface  while  removing  the  bothersome  indentations. 
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5.  DISCUSSION 


It  appears  that  the  general  features  of  an  MHD  solution  for  the  solar  corona  can 
be  simulated  quite  well  by  use  of  an  appropriately  chosen  non-spherical  source  surface. 
Moreover,  the  non-spherical  source  surface  seems  to  generate  a more  satisfactory 
interplanetary  B field  than  do  the  more  conventional  spherical  source  surfaces  of 
Schatten  et  6d.  (1969)  and  Altschuler  and  Newkirk  (1969).  We  take  this  opportunity  to 
discuss  some  of  the  more  subtle  issues  that  may  arise  in  the  course  of  implementing  our 
model  based  on  a non-spherical  source  surface. 

One  might  ask  why  we  have  considered  only  the  internal  sources  of  j^in  specifying 
the  isogauss  source  surface,  when  it  would  have  seemed  more  self-consistent  to  select  an 
isogauss  of  the  ultimate  field  configuration.  We  can  respond  to  this  question  in  several 
ways.  First,  we  note  that  our  results  for  the  dipole  case  are  very  satisfactory,  when 
compared  with  the  MHD  solution  of  Pneuman  and  Kopp  (1971a,b).  Second,  we  remark 
that  the  search  for  a self-consistent  isogauss  source  surface  would  have  entailed  an 
infinite  series  of  iterations  that  would  have  defeated  our  objective  of  defining  an  easily 
implemented  procedure  for  calculating  B.  Third,  we  remark  that  the  minimization  of  tr 
in  (3)  produces  a trajectory  of  extremely  small  (ideally  vanishing)  total  ^ on  the  source 
surface,  i.e.,  at  the  field  cusp  (see  Figures  1-2)  that  constitutes  the  inner  edge  of  the 
neutral  sheet.  Therefore,  the  self-consistent  isogauss  source  surface  (if  definable  at  all) 
is  certainly  not  everywhere  convex,  as  would  be  required  for  the  construction  of  our 
exterior  solution. 


Our  source  surface  bears  a morphological  and  conceptual  resemblance  to  the  self- 
consistent  Alfvenic  surface  calculated  in  MHD  by  Pneuman  and  Kopp  (1971a,b).  Their 
Alfvenic  surface  (at  which  pu  /2  = B /Sir  for  a solar  wind  of  mass  density  p and  bulk 
velocity^),  however,  was  characterized  by  an  equatorial  indentation,  and  so  was  not  a 
convex  surface.  Our  source  surface  rg(0)  is  located  sunward  of  the  MHD  Alfvenic 
surface  r^(  0 ) at  all  latitudes,  i^,  is  completely  enclosed  by  the  Alfvenic  surface  of 
Pneuman  and  Kopp  (1971a, b).  It  is  difficult  to  assess  the  importance  of  the  inequality 
rg(6)<  since,  as  was  stated  above,  we  do  not  know  the  fundamental 

mathematical  relationship  (if  any)  between  source-surface  constructions  and  MHD.  The 
computational  advantage  of  our  source-surface  approach  is  that  we  can  construct  coronal 


and  interplanetary  B fields  from  the  solar  field  without  having  to  calculate  the  solar- 


wind  properties  simultaneously.  However,  we  note  that  Pneuman  and  Kopp  (1971a, b) 


have  defined  r^(e)  without  regard  for  the  thermal  energy  density  contained  in  the 
plasma.  A differently  defined  r^(B)  might  have  had  greater  significance  in  MHD  and 
might  also  have  more  nearly  approximated  our  rg(  0 ). 


One  evident  use  of  the  results  obtained  above  for  the  dipole  case  is  the 
construction  of  a model  solar  wind  in  a prescribed  B-field  geometry.  This  has  been  done 
in  exospheric  theory  (Lemaire  and  Scherer,  1971)  for  a prescribed  radial  magnetic  field 
and  in  fluid  theory  (Durney  and  Pneuman,  1975)  for  a prescribed  non-radial  magnetic 
field,  but  with  solar  rotation  neglected.  It  would  be  highly  relevant  to  develop  model 
solar  winds  in  both  exospheric  theory  and  fluid  theory  for  the  present  model  field,  which 
might  be  viewed  as  a prototype  of  reality.  The  present  model  field  might  also  serve  as  a 
prototype  for  the  description  of  cosmic-ray  transport  within  the  heliosphere.  Jokipii  et 
al.  (1977)  have  considered  the  effects  of  gradient-curvature  and  other  adiabatic  drifts 
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(e.g.,  Northrop,  1963)  on  cosmic-ray  transport  in  a model  field  that,  unlike  (31),  closely 
resembles  the  present  one  at  large  distances  from  the  sun. 


1 


[ 


We  can  make  a qualitative  comparison  between  interplanetary  ^field  observations 
and  our  results  shown  in  Figure  6 by  recognizing  that  our  planar  neutral  sheet 
corresponds  in  nature  to  a neutral  sheet  that  is  warped  about  ±15°  with  respect  to  the 
heliographic  equator  (Schulz,  1973;  Levy,  1976;  Smith  et  al.,  1978).  Thus,  we  should 
interpret  X.  = | ( -n/l)  - 8 | in  Figure  6 as  the  angular  distance  of  the  observer  from  the 
warped  neutral  sheet  at  r ~ 1 AU  and  heliographic  longitude  f , and  we  should  interpret 
\ = 0 as  the  crossing  of  a "sector  boundary"  of  the  sort  described  by  Wilcox  and  Ness 
(1965).  In  a superposed-epoch  analysis  of  their  observations  of  B,  Wilcox  and  Ness  (1965) 
found  a relative  minimum  in  B just  past  the  middle  of  each  sector  (i.e.,  near  the 
maximum  attainable  X ) and  relative  maxima  in  B just  before  and  just  after  the  crossing 
of  a sector  boundary.  A maximum  attainable  \ of  15°  would  correspond  to 
cos  0 w 0.26  in  Figure  6.  Thus,  properly  interpreted,  the  observations  of  B are  in  good 
qualitative  agreement  with  Figure  6.  The  maxima  in  B observed  by  Wilcox  and  Ness 
(1965)  were,  however,  neither  as  pronounced  nor  as  near  the  "sector  boundaries"  nor  as 
symmetrical  about  the  minimum  in  B as  one  would  expect  from  Figure  6.  Thus,  the  non- 
spherical  source  surface  specified  by  (4)  appears  to  account  for  the  qualitative  form  of 
magnetic-field  structure  observed  within  an  interplanetary  sector,  but  the  present 
approach  must  be  applied  to  a more  realistic  model  of  the  solar  B field  before  any 
quantitative  comparisons  can  be  made. 

Similar  superposed-epoch  analyses  reported  by  Wilcox  (1968)  suggest  nearly  a 
doubling  of  the  solar-wind  velocity  (from  u ~ 250  km/sec  to  u ~ 450  km/sec)  between  the 
sector  boundaries  ( X = 0)  and  the  middle  of  a sector  ( X w 15°  by  the  above  inter- 
pretation). This  would  correspond  to  a latitudinal  gradient  - 10-15  km/sec  per  degree  of 
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equivalent  heliomagnetic  latitude  and  is  in  good  agreement  with  the  results  of  Rhodes 
and  Smith  (1976),  who  infer  a gradient  ~ 10-15  km/sec  per  degree  of  heliographic  latitude 
from  spacecraft  observations  made  in  and  near  the  ecliptic  plane.  The  present 
interpretation  suggests  that  numerical  uncertainties  and  scatter  in  the  analysis  of  such 
spacecraft  data  could  be  reduced  greatly  by  taking  account  of  the  observer's  position 
within  a sector,  i^,  of  the  inferred  angular  distance  between  the  observer  and  the 
warped  neutral  sheet. 

In  summary,  we  consider  the  use  of  a non-spherical  (isogaussian)  source  surface  to 
be  very  promising  as  a method  of  modelling  the  coronal  and  interplanetary  B field.  We 
have  demonstrated  the  utility  of  the  method  by  constructing  the  prototype  of  a model 
field,  i^,  that  corresponding  to  an  internal  dipole  surrounded  by  the  source  surface 
specified  by  (4).  We  anticipate  the  application  of  our  method  to  more  realistic 
representations  of  the  solar  B field,  but  we  recognize  that  considerable  work  remains  to 
be  done  before  this  goal  is  fully  achieved. 
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THE  IVAN  A.  GETTING  LABORATORIES 

The  Laboratory  Operations  of  The  Aerospace  Corporation  is  conducting 
experimental  and  theoretical  investigations  necessary  for  the  evaluation  and 
application  of  scientific  advances  to  new  military  concepts  and  systems.  Ver- 
satility and  flexibility  have  been  developed  to  a high  degree  by  the  laboratory 
personnel  in  dealing  with  the  many  problems  encountered  in  the  nation's  rapidly 
developing  space  and  missile  systems.  Expertise  in  the  latest  scientific  devel- 
opments is  vital  to  the  accomplishment  of  tasks  related  to  these  problems.  The 
laboratories  that  contribute  to  this  research  are: 

Aerophysics  Laboratory:  Launch  and  reentry  aerodynamics,  heat  trans- 
fer, reentry  physics,  chemical  kinetics,  structural  mechanics,  flight  dynamics, 
atmospheric  pollution,  and  high-power  gas  lasers. 

Chemistry  and  Physics  Laboratory:  Atmospheric  reactions  and  atmos- 
pheric  optics,  chemical  reactions  in  polluted  atmospheres,  chemical  reactions 
of  excited  species  in  rocket  plumes,  chemical  thermodynamics,  plasma  and 
laser-induced  reactions,  laser  chemistry,  propulsion  chemistry,  space  vacuum 
and  radiation  effects  on  materials,  lubrication  and  surface  phenomena,  photo- 
sensitive materials  and  sensors,  high  precision  laser  ranging,  and  the  appli- 
cation of  physics  and  chemistry  to  problems  of  law  enforcement  and  biomedicine. 

Electronics  Research  Laboratory:  Electromagnetic  theory,  devices,  and 
propagation  phenomena,  including  plasma  electromagnetics;  quantum  electronics, 
lasers,  and  electro-optics;  communication  sciences,  applied  electronics,  semi- 
conducting, superconducting,  and  crystal  device  physics,  optical  and  acoustical 
imaging;  atmospheric  pollution;  millimeter  wave  and  far-infrared  technology. 

Materials  Sciences  Laboratory;  Development  of  new  materials;  metal 
matrix  composites  and  new  forms  of  carbon;  test  and  evaluation  of  graphite 
and  ceramics  in  reentry:  spacecraft  materials  and  electronic  components  in 
nuclear  weapons  environment;  application  of  fracture  mechanics  tc  stress  cor- 
rosion and  fatigue-induced  fractures  in  structural  metals. 

Space  Sciences  Laboratory:  Atmospheric  and  ionospheric  physics,  radia- 
tion from  the  atmosphere,  density  and  composition  of  the  atmosphere,  aurorae 
and  airglow;  magnetospheric  physics,  cosmic  rays,  generation  and  propagation 
of  plasma  waves  in  the  magnetosphere:  solar  physics,  studies  of  solar  nnagnetic 
fields;  space  astronomy,  x-ray  astronomy;  the  effects  of  nuclear  explosions, 
magnetic  storms,  and  solar  activity  on  the  earth's  atmosphere,  ionosphere,  and 
magnetosphere;  the  effects  of  optical,  electromagnetic,  and  particulate  radia- 
tions in  space  on  space  systems. 
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